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Abstract 



We analyze supersymmetric contributions to the branching ratios and CP 
CN '. asymmetries of B^ -^ D^K^ and B^ -^ D^K^ processes. We investigate 

^ . the possibihty that supersymmetric CP violating phases can affect our de- 

^ . termination for the angle 7 in the unitary triangle of Cabibbo-Kobayashi- 

Maskaw mixing matrix. We calculate the gluino and chargino contributions 
^ ■ to 6 ^ u{cs) and h -^ c{us) transitions in a model independent way by using 

O ■ the mass insertion approximation method. We also revise the D^ — D^ mixing 

O ' constraints on the mass insertions between the first and second generations of 

^ ' the up sector. We emphasize that in case of negligible D^ — Z)° mixing, one 

I I should consider simultaneous contributions from more than one mass insertion 

Q \ in order to be able to obtain the CP asymmetries of these processes within 

'^ ' their la experimental range. However, with a large D° — D" mixing, one finds 

a significant deviation between the two asymmetries and it becomes natural 
r^ \ to have them of order the central values of their experimental measurements. 



1 Introduction 

Recently, the BaBar collaborations have measured the charge CP asymmetries Acp^ and 
the branching ratios Rcp± of the B~ -^ D^K^ and B^ -^ D^K^ decays [1]. The following 
results have been reported: 

Acp+ = 0.35 ± 0.13(stat) ± 0.04(syst), Acp_ = -0.06 ± 0.13(stat) ± 0.04(syst), (1) 
R^p^ = 0.90±0.12(stat)±0.04(syst), Rcp_ = 0.86 ± O.lO(stat) ± 0.05(syst). (2) 

These results, with all other 5-factories measurements, provide a stringent test of the 
Standard Model (SM) picture of flavor structure and CP violation and open the possibility 
of probing virtual effect from new physics at low energy. 

In the SM, CP violation arises from complex Yukawa couplings which lead to the 
angles a,l3 and 7 in the unitary triangle of Cabibbo-Kobayashi-Maskawa (CKM) quark 

mixing matrix. The angle /3 = arg ( — y'^y^ ) has been determined by the CP asymmetry 
in B^ -^ J/ipKs process which is dominated by tree level contribution. Concerning 
the angle 7 = arg ( — y"'^y"'' ), it is believed that a theoretically clean measurement of 
this angle can be obtained from exploiting the interference between B^ -^ D^K^ and 
B^ -^ D^K~ when D^ and Z)° mesons decay to the same CP eigenstate [2]. 

At the quark level, the B^ -^ D^K^ and B^ -^ D^K^ decays are based on 6 ^ u{cs) 
and h -^ c{us) transitions respectively. Therefore, their SM contributions at tree leve are 
suppressed by the CKM factors VcsV*^, and V*gVch which are of order 10"^. This gives the 
hope that it may be possible for a new physics beyond the SM, like supersymmetry, which 
contributes to these decays at one loop level to manifest itself and compete the SM. In 
this paper we aim to investigate this possibility and check, in a model independent way, 
whether supersymmetry can significantly modify the CP asymmetries in B~ — > DK~ 
processes and hence affects the determination of the angle 7. Therefore, we perform 
a systematic analysis of the SUSY contributions to B ^ DK processes. We compute 
SUSY contributions to 6 ^ u{cs) and h -^ c{us) transitions through the gluino and 
chargino exchange, using the mass insertion approximation method. This approximation 
is quite useful tool for studying the SUSY contributions to the flavor processes in a 
model independent way. We show that the gluino box diagrams give the dominant SUSY 
contribution while the chargino exchanges lead to subdominant contributions. 

It turns out that the D^ — D^ mixing may limit the gluino contribution to B^ -^ DK~ 
due to the stringent constraints on the mass insertions between the first and second 
generations in the up sector, (5^s)i2- Thus in our analysis, we revise the D^ — D^ mixing 
constraints [3] and take them into account. We find that with a single mass insertion, 
the SUSY contribution to B^ -^ DK~ decay will be much smaller than the SM result. 
Nevertheless, with simultaneous contributions from more than one mass insertion, the 
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SUSY effect can be enhanced and the results of the CP asymmetries become within la 
experimental range, while the D^ — D^ mixing constraints are satisfied. 

The paper is organized as follows. In section 2 we study the CP asymmetries and the 
branching ratios of B~ -^ DK~ in the SM. We show that in the SM the branching ratios 
Rcp± are within the experimental range. While the CP asymmetry Acp+ is below its la 
experimental lower bound and the value of Acp_ is typically negative. In section 3 we 
compute the gluino and chargino contributions to 6 ^ u{cs) and h — * c{us) transitions in 
terms of the mass insertions. Section 4 is devoted for analyzing the SUSY contribution 
to D^ — Z)° mixing and revise the possible constrain on the mass insertions (5^^)i2- 
The analysis of SUSY contribution to CP asymmetries Acp^ and branching ratios Rcp± 
is given in section 5. We show that in case of negligible D^ — D^ mixing, one should 
consider simultaneous contributions from more than one mass insertion in order to obtain 
Acp± within their la experimental range. Nevertheless, the usual relation: A+ ~ — A_ 
which is valid in the SM remains hold. With a large D^ — D^ mixing, one finds a significant 
deviation between A^ and A_ and it becomes natural to obtain Acp^ of order the central 
values of their experimental measurements. Finally, we give our conclusions in section 6. 

2 B ^ DK in the Standard Model 

In this section we analyze the CP violation in B~ -^ DK^ decays within the SM. The 
possible quark level topologies of B^ — > DK'^ that contribute to the amplitude A{B^ -^ 
D^K~) and A{B~ -^ D^K~) in the SM can be classified to the following three categories, 
as shown in Fig.(l): color- favored tree (T), color- suppressed tree (C) and annihilation 
(A). These processes are given in terms of the CKM factors Ac = VchV*^, \u = VubV*^. 
The decay B^ -^ D^K~ receives contributions from T and C with factor Ac, while 
B^ -^ D^K^ get contributions from C and A in terms of A^. Since the contributions 
from the annihilation process to the matrix elements are quite suppressed at the leading 
order correction [4], it is quite reasonable to assume that A = 0. In our analysis we will 
adopt this approximation and therefore the general parametrization of the SM amplitudes 
of B^ -^ DK^ decays can be given by 

A^^{B- ^D^K') = \Ai\e'^-' =f + C, (3) 

A^^{B- ^D^K-) = \A2\e'^^e'^ = C, (4) 

where 6i, 2 = 1,2 are the strong (CP-conserving) phases. T and C refer to the color 
allowed and color suppressed tree amplitudes involving b -^ c{us) while C is related to 
the process b -^ u{cs). In terms of the two CP-eigenstates of the neutral D meson system, 
D^jp^ = (D^ ± Z)°)/v^, one considers the ratios Rcp± of charged averaged partial rates 
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Figure 1: SM contributions to i? -^ DK : color- favored tree (left up), color-suppressed 
tree (right-up and left-down) and annihilation (right-down) 



and the charge asymmetries Acp^' 
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We define the ratio of the SM amplitudes of i? -^ D^K and B -^ D^K as 

AS^{B- -^ D^'K-) 



rpe'^^e'^ 



ASM^B- -^ D^K- 



(7) 



According to Eqs.(3,4), r^ = IA2/A1I and Sp = S2 — Si. Using this parametrization, one 
finds that R± = Rcp± is given by 



R± = 1 + Tp ± 2rB cos 6b cos 7, 



and A± = Aqp^ takes the form 

±2rB sin dp sin 7 



A, 



From Eq.(8) one gets 



±2rB sin Sp sin 7 



1 + r^ ± 2rB cos 5p cos 7 



i?, 



CP± 



cos 7 



R^ — R— 

Atb cos 5b 



(8) 



(9) 



(10) 



Thus, by using the expressions for the CP asymmetries A± in Eq.(9), one can factorize 
the dependence on the strong phase and gets the following expression for the angle 7 in 
terms of R±, A± and tb only: 

2cos7(A+-A_) R+R_ 
sm7 = = . (11) 

'l6r|cos2 7-(i?+-i?_)2^+ + R- 



From this expression, one can easily see that the central experimental values of R± and 
A± with Tb — 0.1 implies that the angle 7 is of order 7 ~ 71°. It is worth mentioning 
that within the SM, the effect of the D^ — D^ mixing is very small on extracting the angle 
7 using the B^ -^ DK^ decays. As emphasized in Ref.[5], neglecting this mixing implies 
an error in determining 7 of order 0.1 — 1°. 

In order to analyze the SM predictions for the A± and R± and compare them with 
the experimental results reported in Eqs.(l,2), let us consider the SM contributions to the 
h -^ u{cs) and h -^ c{us) transitions. As shown in Fig. 1, within the SM the B^ -^ DK~ 
are pure 'tree' decays. The effective Hamiltonian of this transition is given by 

Hlr^''^ = ^Kbv;: [cMQi + c^2(/i)g^] , (12) 

where Ci and Q^ are the Wilson coefficients and operators of this transition renormalized 
at the scale fj, with 

g^ = (m"7^^&°) {s^rLcf") , Q^ = (u'^i^^Lb^) (A^^c") , (13) 

where L = (1 — 75). The effective Hamiltonian for the b -^ c{us) transition can be 
obtained from the effective Hamiltonian in Eq.(12) by exchanging u ^^ c. The SM results 
for the corresponding Wilson coefficients are: 

C,{mw) = 1 - ^^, C,{mw) = ^. (14) 

However, due to the QCD renormalization to the scale /i ~ m?,, Ci and C2 get mixed, as 
will be discussed in more details in the next section, and one finds 

Ci(/i) = 1.07, C2{^l) = -0.17. (15) 

To evaluate the SM results to the decay amplitude of B^ -^ DK^ , we have to de- 
termine the matrix elements for the operators Q" 2- A detailed analysis for the matrix 
elements will be given in the next section. Here, we just give the matrix elements for 
these four operators in naive factorization: 

{D'K-IQ'IIB-) = -^, {D'K-mB-) = -X, (16) 
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{D'^K-lQllB-) = -\x-Y, (D°ir-|g^|5-) = -X-iF, (17) 



where 



B^D, 



X = iF^~-^{m'j,)fn{m'B ' O, ^ = zFo^^^(m^)/^(m^ - mj, 



(18) 



There are two comments in order: i) The naive factorization can not determine the 
strong phases, therefore, in our analysis we consider these phases as free parameters, 
ii) As mentioned above, the factorized matrix element {D^K~\{s'y'^Lc)\0){0\u'jfj,Lb\B~) 
corresponding to an annihilation process is suppressed as showed in Ref. [4], and can be 
neglected. Therefore, 
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3 / V 3 

Fixing the hadronic parameters as follows: /d = 0.2, fx = 0.16, Fq^^ = 0.34, Fr 
0.62, and the meson masses as: mx = 0.49, rriD = 1.86, and rriB 

tb ^ 0.05 



(19) 
(20) 



B^K 



5.278 GeV. One finds 



(21) 



Note that it is customary assumed that with a large uncertainty, the SM prediction for 
rB may be much larger than the above value (can be (9(0.1), see Ref.[6]). Here we will 
use the value that we obtained in Eq.(21) as a typical value for the SM contribution. In 
order to have a general picture of the SM predictions for the CP asymmetries A± and the 
branching ratios R±, we plot in Fig. 2 R± versus A±. Here, we vary the parameter 6b in 
the range [0, vr] and the angle 7 is also considered between and n. 
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Figure 2: i?+ versus y4+ and R- versus A^ within the Standard Model. The horizontal 
line in the left figure represents the lower bound of A^ at la experimental range. 



As can be seen from the results in Fig. 2, the SM predictions for the branching ratios 
R± are within the la experimental range. However, the results for the CP asymmetry A+ 



are below its experimental lower bound. Also the SM leads to a negative CP asymmetry 
A_ which is still consistent with its experimental results in Eq(l), due to the large un- 
certainties in these measurement. Therefore, more precise measurements would be very 
important in analyzing the SM predictions for Rj- and A± and, hence, in determining the 
value of the angle 7. 

3 SUSY contributions to 6 ^ u{cs) and b -^ c{us) 

The crucial point to note from the previous section, is that the SM contributions to the 
amplitudes oi b ^ u{cs) and b -^ c{us) transitions are suppressed by the CKM factors 
Vub — 0(10^^) and V*gVcb — 0(10^^) respectively. Therefore, it may be possible to have 
a comparable effect from new physics at one loop level which can compete with the SM 
tree level contribution. In this section we study the supersymmetric contributions to the 
b -^ u{cs) and b -^ c{cs) transitions. In this case, the effective Hamiltonian H^^^ for 
the b -^ u{cs) can be expressed as 

10 

^cff''=' = E {cr{^^) gr(/i) + An-") Qr(/i)) , (22) 

where C" are the Wilson coefficients and Q'^ are the relevant local operators at low energy 
scale IX ~ iTLiy. The operators Q" are given by 



Ql = (m"L6") (s^Lc'^j , Ql = (vf'Lb'^) [s'^Lc'') , (23) 

g^ = (m"L6") (sf^Rc^) , Ql = (u'^Lb^) (sf^Rc" 



Ql = (txV„,L6") (s'^a^'^Lc^) , g^o = (izV„,L6'3) (s'^a'^'^Lc 



where a and (3 refer to the color indices. L,R are given by (1 =i= 75) respectively and 
a^'^ = ^[7'^, 7^^]. The operators g" are obtained from Q^ by the chirality exchange L ^^ R. 
In the SM, the coefficients Cf are identically vanish, while in SUSY models, they receive 
contributions from both gluino and chargino exchanges. The corresponding operators for 
b -^ c{us) can be obtained from the above expression by exchanging u *^ c. 

The dominant SUSY contribution to the b -^ u{cs) transition can be generated through 
the box-diagrams with gluino exchange, as in Fig. 3, and chargino exchange, as in Fig. 4. 
From these figures, one can see that the b — > u{cs) transition is based on two topologically 
distinct box diagrams only for gluino or chargino exchange. This is unlike the b -^ d and 
6 — s> s transitions that contribute to B — B mixing, where four topologically distinct box 
diagrams are included [7]. Therefore, it is expected that the Wilson coefficients for this 
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Figure 3: Box diagrams for i? ^ K D^ {b -^ u{cs) transition) with gluino exchanges, 
where h, k,m,n = {L, R}. 
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Figure 4: Box diagrams for B ^ K D^ {b -^ u{cs) transition) with chargino exchanges, 
where U = {u, c, t}, D = {d, s, b} and h, k,m,n = {L, R}. 



process are different from those obtained in the hterature for 6 ^ s transition. It is 
also worth mentioning that contributions through penguin diagrams to these transitions 
are always hybrid {i.e., contain internal SUSY and SM particles). Therefore, they are 
suppressed by Vub in addition to the usual loop suppression factor, hence they are much 
smaller than the pure SM or pure SUSY contributions. Thus, the Wilson coefficients at 
mw scale can be expressed as follows 



c- = (cf)*'^ + (ci^y + (cf )^ 



(24) 



We evaluate the SUSY contributions to the Wilson coefficients by using the mass 
insertion approximation. The Mass insertion approximation is quite useful method in 
order to perform model independent analysis of flavor changing processes in general SUSY 
models. In our analysis we set to zero the contributions that are proportional to to the 
Yukawa coupling of light quarks. Also, we use the approximation of retaining only terms 
proportional to order A. In the case of the gluino exchange all the above operators give 
significant contributions and the corresponding Wilson coefficients are given by 
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Cl{mw) 






-^MlMSIl)^ [Mx] 



Axfeix) 
20a; /e (a; 
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(26) 



'±oTfl 



Cl{mw) = 



Climw) = ^(4L)23(5k)i2:r/6(x), 



'limw) = ^{SL)23{S'kL)i2xMx), 






^,{Kr)23{S1l)i2[-Mx) + 7xf,{x) 



12m 
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Climw) = ^,{Kr)23{51l)i2[5Mx)+xUx) 

4o?7i 
C~'oi^w) = -^Kl)23Kl)i2xMx)- 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
(34) 



where x = rn^/fh^. The nig is the gluino mass and the m^ is an average squark mass. The 
functions f%{x) and f%{x) are the same as the loop function obtained in case of 6 — *> d{qq) 
and are given by 



_ 6(l + 3x)lnx + x3 -9x2-9x+ 17 
6(a; — 1)^ 
~ 6a;(l + x) Inx — a;^ — 9x^ + 9a; + 1 
^^^") = 3(^^I? • 



(35) 
(36) 



The Wilson coefficients Cf are simply obtained by interchanging L ^^ J? in the mass 
insertions appearing in Cf . The above Wilson coefficients are due to the gluino exchange 
of 6 -^ -u transition, the corresponding coefficients for b ^ c transition can be obtained 
by changing the mass insertions (5ab)i2 fo (^^5)21 where {A,B} = {L,R}. 

Note that the discrepancy between the above SUSY Wilson coefficients of 6 ^ -u 
transition and those of 6 ^ rf or 6 — * s AB = 2 transition is due to the following reasons: 
1) the b ^ u transition is based, as mentioned above, on two distinct box diagrams only 
in contrast of the AB = 2 transition where four distinct box diagrams are involved. 2) All 
the external quarks in the box diagrams of 6 -^ u{cs) are different, therefore, one can not 
use Fierz transformation to relate any operator with the other unlike the case in AB = 2. 
For instance, in Bd — Bd mixing the operate Q2 = {d°''y^Lb^] (d^'jf^Lb"] is equivalent to 
the operator Q2 = (d^'^^Lb") (d^^y^Lb^). In this case, the Wilson coefficients Ci and C2 
in Eqs. (25,26) are combined together and leads to to the usual AB = 2 Wilson coefficient: 
Ci oc as/{l08m'^) f24a;/g(x) + 66/6(a;) j [7]. In this respect, it is clear that the expression 

used in Eq.(9) in Ref. [3] for H^^ is incorrect. 



Now let us turn to the chargino contributions to the effective Hamiltonian in Eq.(22) 
in the mass insertion approximation. The leading diagrams are illustrated in Fig. 4, where 
the cross in the middle of the squark propagator represents a single mass insertion. Within 
the above mentioned approximation where we neglect contributions proportional to the 
light quark masses, one find that the relevant chargino exchange affects only the operator 
Qi, as in the SM and the corresponding Wilson coefficient is given by 



CU^w) 



a\ a r 
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V^v;,v,,UaU;,{6iL)u mL)2z + Ks 



LLJIS) 



'An 



U^iUnV,rV:^{bi^),^ m^hs + KSl^)n) 



{L2{xi,Xj) - 2Lo{xi,Xj)) ,(37) 



where a = g'^/An and g is the SU{2) gauge coupling constant. The A parameter stands 
for the Cabibbo mixing, i.e., A = 0.22. The Uij and Vij are the unitary matrices that 



;^2 



diagonalise the chargino mass matrix and yt is the top yukawa coupling. The Xi = rn^Jm 
and the functions Lo{x,y) and L2{x,y) are given by [7] 



Lo{x,y) = 


= ^/xy 


x-y 


L2{x,y) = 


xh2{x) - 


-yh2{y) 



x-y 



ho{x) = - 
h2{x) 



{l-xf 

2 + 3a; - a;2 



+ 



X) 



6 In a; 

{1-xy 
6a; In a; 

d-xY' 



(38) 



Finally, we have also neglected the small contributions from the box diagrams where both 
gluino and chargino are exchanged as in Ref . [8] . 

To obtain the Wilson-coefficients at the scale rrib one has to solve the corresponding 
renormalisation group equations. The solution is generally expressed as 



Ciirub) =J2Uij{mb,mw)Cj{mw), 



(39) 



where Uij{mb, mw) is the evolution matrix given by the 8x8 anomalous dimension matrix 
of leading order (LO) corrections in QCD [9]. Note that we have not included the operators 
Qg^io since they have zero matrix elements at LO and also they do not mix with the other 
operators in the evolution from mw scale down to rrib scale. 



/. 



U{mb,mw) = V 
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-1 



(40) 
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D 



where V diagonalizes the 7^^)^ 

and 7*^°) is the diagonal elements of Xd ■ The value of (3o is given by (3o = yA^'c — |/ 
where Nc is the number of colors and / is number of active flavors. Finally, the anomalous 
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dimension matrix 7° at the leading order is given by 
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(42) 



As can be seen from the above matrix that the mixing between different operators is 
divided into blocks. Each block contains two operators {Qi,Qi+i), i = 1,3,4,7 and with 
no mixing between different blocks [10]. 

Let us now consider the evaluation of the hadronic matrix elements of the above 
operators which represents the most uncertain part in this calculation. In the limit of 
neglecting QCD corrections and rrib ^ ^qcd, the hadronic matrix elements of B^ -^ 
DK~ decay can be factorized. The hadronic matrix elements for the operators Q" are 
given by 



(43) 



and 
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"e" 

X 

y 

{D'K-\Q'i,\B-)=Q. 



(44) 



While the hadronic matrix elements for the operators Q1 are given as follows: 

{D'K-\Q'i\B-) = Q, 



(45) 



and 



{D''K-\ Ql \B-) 



3 ' 
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(46) 



where X and Y are given in Eq.(18). 

Having evaluated the SUSY contributions to the Wilson coefficients and determined 
the hadronic matrix elements of the relevant operators, one can analyze the decay ampli- 
tude of B~ -^ DK~ and study the SUSY effect on the CP asymmetries A± and branching 
ratios R±. As can be observed, the Wilson coefficients depends on several mass insertions 
which are in general complex and provide new sources for the CP violation beyond the 
SM phase in the CKM mixing matrix. These new CP violating phases may contribute 
significantly to the b -^ u transition and affect the determination of the angle 7. Nev- 
ertheless, one should be very careful with the constraints imposed on these parameters. 
In fact, the dominant gluino contributions depend on the mass insertions: ((5^b)i2 and 
(5^^)23. The mass insertions ((5^^)23 are constrained by the experimental results for the 
branching ratio of -B — > Xs'~f [11-13]. These constraints are very weak on the LL or RR 
mass insertion and more stronger for LR or RL mass insertions. Concerning the mass 
insertion ((5^^) 12, the chargino contributions to the K^ — K'^ impose constraint on the 
LL mass insertion only [14]. However, the D^ — Z)° mixing may induce more strangest 
constraints on both LL{RR) and LR{RL) mass insertions. Therefore, before we proceed 
in analyzing the decay amplitude of B^ -^ DK^ , we will take a short detour and give a 
detail analysis for the SUSY contributions to D^ — D^ mixing and revise the corresponding 
constraints on the {6"^^) 12 mass insertions. 

4 Constraints from D^ — D^ mixing 

We start this section by summarizing the SM results for the D^ — Z)° mixing, then we con- 
sider the supersymmetric contribution to the effective Hamiltonian for AC = 2 transitions 
given by gluino and chargino box exchanges. 

In the D^ and D^ systems, the flavor eigenstates are given by D° = (uc) and Z)° = (uc) . 



11 



The formulism for D^ — D^ mixing is the same as the one used for K^ — K^ and B^ — B^ 
mixing. The mass eigenstates are given in terms of the weak eigenstates as: 

Di,2 =pD^±qD\ (47) 

where the ratio q/p can be written in terms of the off-diagonal element of the mass matrix: 



?/P = y MI2/M12 and q/p 7^ 1 is an indication for the CP violation through mixing. The 
strength of D^ — D^ mixing is described by the mass difference 

^Md = Md,-Md,. 

The present experimental results for AMd is given by [15] 

(AMz5),,p < 1.7 X 10-^-^ GeV. (48) 

The CP asymmetry of the D^ and Z)° meson decay to CP eigenstate / is given by 

r(DO ^ /) - r(DO ^ /) 



a fit) 



Sf sin(AMDt) + Cf cos(AMz)t), (49) 



2Im 

Sf 



where Sf and Cf represent the mixing and direct CP asymmetry respectively and they 
are given by 

|p(/)P + i' ^' - M7)FTT- ^'°^ 

The parameter p(/) is defined by p(/) = Jrj^o^fl ■ Generically, the AMc and Sf can be 
calculated by 



AM 



D 



{D'\H^fi='\D') , Sf = sin (arg {D'\H^/'='\D') . (51) 



Here H^f'j ^ is the effective Hamiltonian responsible for AC = 2 transition. In the 
framework of the SM, this transition occurs via box diagram in which two virtual down 
quarks and two virtual W bosons are exchanged. The AM|^ = 2|(L'°|(ifg^J=2)'5^|L)0)| 
is given by [16] 

AM^^ =. ^,Mlfl\V:,V^\' ^'^' ~ f'^' ^ 1.4 X 10-i« GeV. (52) 



37r2" 



m 



As can be seen from this expression, the SM predicts a very small D^ — Z)° mixing. Note 
that, in the above estimation for AM^^ , the 6-quark contribution has been neglected 
since it is much smaller due to the CKM suppression. Also, the CP violation is absent 
in the mixing and in the dominant tree level decay due to the involving of the first two 
generations only. 
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In supersymmetric theories, the dominant contributions to the off diagonal entry in 
the D^-meson mass matrix, 7Wi2 = {D^\H^^'^\D^) , is given by 



M 



12 



M 



SM 
12 



M' 



12 



M^ 



12 5 



(53) 



where M.12-, and A^i2 correspond to the gluino and chargino contributions respectively. 
The effect of supersymmetry can be parameterized as follows 



^2g2jec 



M 



12 



where AM^ 



2|A1^M|^2 ^^^ 29^ 



arg 1 



Mf^' 



A^fusY 



Mff 



(54) 



As in the case of K^ and B^ 



systems. The most general effective Hamiltonian for AC = 2 processes, induced by gluino 
and chargino exchanges through box diagrams, can be expressed as 



H. 



AC=2 
eflf 



J2 Ci(/i)Qi(/i) + E Ciifi)Qiifi) + h.c. , 



(55) 



j=i 



i=l 



where Ci{fj,), Ci{fi) and Qi{fi), Qi{lj) are the Wilson coefficients and operators respectively 
renormalized at the scale /i, with 



Qi 

Q4 



-a a -13 P 

-a a -13 13 



Q2 

Q5-- 



-a a -13 13 
"^R^L "^R^Li 

-a P -13 a 

^r'^l "^l^r ■ 



Qa 



-a 13 -13 a 
^R'^L "^R^Li 



(56) 



In addition, the operators <5i,2,3 are obtained from <5i,2,3 by exchanging L *-^ R. 

In the case of the gluino exchange all the above operators give significant contributions 
and the corresponding Wilson coefficients are given by [17] 



Ci{mw) 
Climw) 
Ci{mw) 



at 



216m? 



at 



216m? 



(24xUx) + 66h{x)){6 



20AxMxM,)l^, 



u \2 
12)LL^ 



at, 



216m? 



36a;/6(a;)(5r2)k, 



(57) 



at 



216m? 



at 



216m? 



(504x/6(x) - 72/6(x)) {5^^2)ll{5'12)rr - l^2Ux){5^^^)LR{5'i^)RL 
(24x/6(x) + 120/6(x)) {5\2)ll{5'^2)rr - im,{xWi2) lr{5'^2)rl 



where x = rnym and m is an average squark mass. The functions f&{x) and Jq^x) are 
given in Eqs. (35,36). The Wilson coefficients 6*1-3 are simply obtained by interchanging 
L *^ Rm the mass insertions appearing in Ci_3. 
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In the case of the chargino exchange the operator Qi only gives a significant contri- 
bution [14]. At the first order in the mass insertion approximation, the Wilson coefficient 
Ci{mw) is give by 



Cf{mw) 



9^ 



7687i'^m 



Zi rt-y-t Zi 



E |V^.inv;-in'^2i)L^2(:r„a:,). 



«j 



where Xi = 7n-+/m , and the function L2{x,y) is as given in Eq.(38). 



X- 



As usual, the Wilson coefficients Ci{j2) are related to Ci{mw) by [18] 



(58) 



(59) 



where rj = as{mw)/cisif^) and the coefficients b^ , c^ , and Oj appearing in (59) can be 
found in Ref.[18]. Also the matrix elements of the operators Qi in the vacuum insertion 
approximation are given by [17] 



{D"\Qi 

{D'\Q2 



{D"\Qt 



|D°) 



,;Mnfn, 



5 



24 \m 



M 



D 



m,, 



MdI 



Di 



1 



24 Vm, 



mM 



D 



m„ 



Mofl, 



(60) 



1 1 
24 4 \m 



M 



D 



m,, 



1 1 



M, 



D 



12 Xrric + rriu 



MdH, 



The same results are also valid for the corresponding operators Qi since strong interactions 
preserve parity. 

We now discuss the results of SUSY contribution to D^ — D^ mixing. It is worth 
mentioning that the mass insertions (5^^)12 are strongly constrained by the experimental 
limits of K^ — K^ mixing. In particular, the AMk upper bound implies that |((5ii,)i2l < 
10~^ [17]. Therefore, the chargino contribution to AM^y becomes very suppressed and can 
be neglected with respect to the gluino contributions which depend on less constrained 
mass insertions (^^5)12 • As an example, we present the gluino contribution to AM^), 
with rrig ~ rrig ~ 500 GeV: 



1.7 X 10-13 



33.4(5 



u \2 
LL) 12 



1733.6(52^)^2 - 3178.5((52^)i2(5]^Ji2 + 1733.6(5]^^ 



u \2 
12 



12M6.9i5l^),2{5''nn)i2 + S3A{5l^ 



u \2 
12 



< 1. 



(61) 
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X 


V (^LL(Ri?))l2 


V (^LR(i?L)) 12 


V {nL)l2{5lR)l2\ 


^J {51r)i2{51l)i2\ 


1/4 
1 
4 


7.5 X 10-2 

1.7 X 10-2 

0.4 


2.1 X 10-2 
2.4 X 10-2 
3.3 X 10-2 


7.5 X 10-3 
8.7 X 10-3 
1.2 X 10-2 


1 X 10-2 

1.7 X 10-2 
4 X 10-2 



Table 1: Upper bounds on 

1/4,1,4. 



m 



AB))12\ 



from AMz) < 1.7 x 10 ^3 GeV, for x = myntj. 



From this expression, we can see that the strongest constraint will be imposed on the 
product (5ll)i2(<^^r)i2 while the constraint obtained on the individual mass insertion 
{51l)i2 or {5'^r)i2 is less stringent. 

As usual in this kind of analysis, the most conservative constraints on the mass inser- 
tions can be obtained by considering the contribution due to a single mass insertion per 
time and set all other ones to zero. In table I, we present the results for the upper bounds 
on the relevant mass insertions from the experimental constraint on AM^ for x = 1/4, 1, 
and 4. We find that these bounds on (5^b)i2 are more stringent than those obtained from 
the chargino contribution to the K^ — K^ in Ref.[14]. In fact, the ((5]^r)i2 and (5rj^)i2 are 
completely unconstrained by the chargino contribution to K^ — K^ mixing. Therefore, 
their bounds in the above table are the only known constraints. However, we should 
mention that these constraints may be relaxed if one consider simultaneous contributions 
from more than one mass insertion. In this case, a possible cancellation may occur which 
reduce the SUSY contribution significantly and leave a room for a larger mass insertion. 

Finally, we comment on the CP violation in this process. As emphasized above, the 
SM contribution to D^ — D^ is real since it is proportional to V*gVcd- Furthermore, it is 
much smaller than the dominant gluino contribution. Therefore CP violating phase 6^ in 
Eq.(54) can be written as 



9, 



I 
c - - arg 



M 



12 



Mff 



- arg [M.^ 



12 



(62) 



In case {Sijiju gives a dominant contribution to A^i2, Oc will be given by 



e, = -e.Tg{{6l^U2)^0{l), 



(63) 



which means that the mixing CP asymmetry of D^ — Z)°, Sf, could be quite large. There- 
fore, one can conclude that the new physics in general and supersymmetry in particular 
could enhance the D^ — D^ mixing significantly. 
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5 Supersymmetric contribution to R± and A± 

In this section we study the supersymmetric contributions to the CP asymmetries and 
the branching ratios of B^ -^ DK~ decay in the following cases: 1) negligible D^ — D^ 
mixing. 2) Large D^ — D^ mixing due to a possible significant SUSY contribution as 
advocated in the previous section. 

In general, applying the naive factorization approximation implies that the amplitudes 
A{B^ -^ DK^) are given by 



A{B- ^ D^K-) = Y: {c^ - C^) {D^K'mB-), 



i=l 



and 



A{B- ^ D'K~) = E (^r - Cr) {D'K-\Q^^\B-). 



(64) 



(65) 



1=1 



The sign difference between the Wilson coefficients Cj and C*j in the above equations is 
due to the fact that the initial and final states of B~ -^ DK~ decays have opposite parity 
and therefore {DK~\ Qi \B-) = -{DK-\ Qi \B-) [19]. 



5.1 R± and A± with negligible D^ — D^ mixing 

In case of neglecting the effect of D^ — D^ mixing, it is useful to parameterize the SUSY 
contribution by introducing the ratio of the SM and SUSY amplitudes as follows: 



and 






Ri e 



i{<t>i~l) ^i-h 



i(p2 „i&2 



Ro e'^^ e''\ 



(66) 



(67) 



ASM^B- -^ D^K-) 

where Ri stands for the corresponding absolute value of \A^^^'^ /A^^\, the angles 0j are 
the corresponding SUSY CP violating phase, and 5i = 5f^ — 5f^^^ are the strong phases. 
In this respect, our previous definition for the SM ratio of the amplitudes of B~ -^ D^K~ 
and B~ -^ D^K~ in Eq.(7) will be generalized as follows 



A{B- 
A{B- 



oii--) 



D^K 



ASM^B- -^ D^K-) + A^usY^B- -^ D^R- 



where 



R 



B 



Tb 



TbC 



e^T + Rie' 



iSf 



t(pi 



1 + i?2e*'^2 



e'^ + Rie 
1 + i?2e^<^2 

and (pB = 
16 



R^e'^B^i<PB^ 



arg 



e'^ + Rie^<^ 

TTr^ 



(68) 



(69) 



Note that, for simplicity, we have assumed that the SM and SUSY strong phases are 
equal. In this case, the ratios R± and the CP asymmetries A± take the form 

R± = l + Rl± 2Rb cos 6 b cos 0b, (70) 

_ ±2Rb sin Sb sin 0b 

1 + R%± 2Rb cos 5b cos 0b ' 
As shown in Eq.(69), the deviation of Rb from the standard model value tb is governed 
by the size of Ri and i?2- Therefore, we start our analysis by discussing the dominant 
gluino contributions to Ri and R2. We choose the input parameters as m = 250 GeV, 
a; = 1 we obtain 

i?i = 0.15((5^J23('5ll)i2-0.17(O23(<5]^Ji2 + 0.18(4J23(5lb)i2-{^ -- ^}, (72) 

and 

i?2 = -0.01(5^^23(5^)21- 0.015(4^)23(52^21 +0.03(4j23('5L^)21-{i^-i?}. (73) 

Using the fact that the mass insertion is less than or equal one, we find that R2 '^ 1, 
i.e A'^^^^iB" -^ D^K~) < A^^^{B- -^ D^K). It is worth mentioning that (51b)23 
is constrained by the experimental results for B -^ Xg'-j decay. These constraints are 
very weak on the LL and RR mass insertions and they can be of order one. How- 
ever, they impose stringent upper bounds on the LR and RL mass insertions, namely 
l(^Li?(RL))23| ~ 1-6 X 10~^ [11]. Concerning the (5^b)i2; the important constraints on 
these mass insertions are due to the D^ — D^ mixing. Applying these constraints one 
finds that Ri is also quite small and the SM gives the dominant contribution. Therefore, 
there will be no chance to modify the results obtained in Fig. 2. However, as advocated 
in the previous section, these constraints can be relaxed if one allows for simultaneous 
contributions from more than one mass insertion, which is the case in any realistic model. 
In this case, there may be cancellation between different contributions which reduces 
the SUSY contribution to D^ — D^ mixing without severely constraining the mass inser- 
tion. If we adopt this scenario and assume, for instance, that (5^^)23 — -(5^^)23 and 
(5ii)i2 — ~('^fi;R)i2, then one can easily see that Ri ~ (9(0.6) and the phase 0i is given 

byarg[(5ij23 + ('5EL)i2]. 

In this case, one can easily observe that different combinations of (7, 0i) can lead 
to values for the A± within the experimental range. Therefore, the supersymmetric CP 
violating phases may affect the extraction of the angle 7. As an example, let us consider 
the case where Rb is enhanced from 0.05 (SM value) to 0.1 and the phase 0b is given by 
70°, which can be obtained by 7 ~ tt/S and 0i ~ 7r/2 or 7 ~ 7r/2 and 0i ~ 7r/3. In this 
case, one finds that 

R+ ~ 1.1, R- ~ 0.94, A+ ~ -A_ ~ 0.2. (74) 
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Therefore, we can conclude that the SUSY contributions to B~ -^ DK~ imply that A+ 
and A_ are within their la experimental range simultaneously, unlike the SM results. 

Finally, it is important to mention that in this scenario it is a challenge to find a 
realistic SUSY model that accommodates these results and satisfies all other constraints. 
Also the observation of A+ indicates that the ratio of the amplitudes for the processes 
B^ -^ D^K" and B^ -^ D^K" is larger than 0.1 which is rather difficult to obtain in 
the SM, so it may be a hint for a new physics effect. 



5.2 R± and A± with large D — D mixing 

In the previous analysis, we have ignored the effect of the D — D^ mixing. Now we consider 
this effect and define the time dependent meson state \Di) = \D^(t)) and ID2) = \D^{t)) 
as 



1^1) = g^it)\D') + ^ g-itW) 

P 

ID2) = g^it)\D') + ^g4t)\D') 



(75) 
(76) 



where q/p is defined, as in the previous section, by 



P 



IM*, 



M 



12 _ g-2jee 



12 



(77) 



As shown in Eq.(63), the phase O^ is of order one . The functions g±{t) is given by [20] 



9± = -^\e 



-Hit 



±e 



-4/12* 



9+{t) 
9 -it) 



^i~iMot-T/2) 



^{-iMot-r/i) 



{-ixD - yDfr/2 + . 



Here r = Ft. In this case, the decay amplitudes of -B -^ DK are given by 

A{B- -^ DiK') = A{B- -^ D°K-)g+{t) + A{B- -^ D^R-) - g^{t), 



and 



A{B- -^ D2K') = A{B- -^ D^K-)g+{t) + A{B- -^ D^K') - g_{t). 



(78) 



with /ij = Mui — iT£,-/2. In terms of Xn = ^^n. g^j^^j y^ _ ^^ where T = T^^ + r2, one 
finds 



Also the decay rates are defined as [20] 

T{B~ ^ DK-) = I dt 
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A{B- 



DR- 



(79) 
(80) 

(81) 

(82) 

(83) 



Therefore, one finds that 

r{B- -^ DiK-) = \A{B -^ D^K-)\'^ (^G+ + RlG^ + 2RBRe [G'+_e-^('^«+'^«-2^^) 

where Gi are given by 



(84) 



oo 



G+ = I \9At)\'dtc.^(^l + yl±^y (85) 

G- = r \9-it)\'dt -^ f [^^^) , (86) 

G^- = r94t)9-m^l[ -y^-'''^ ). (87) 



The CP asymmetries Acp^^ S'-re defined by 

, ^ nB- -^ Di,2K-) - r(i?+ ^ Di,2ir+) 

^^1.2 r(5-^Di,2ir-) + r(5+^Di,2i^+)' ^ ^ 

Thus one can easily prove that 

, Rb [yp sin 6b sm{(f)B - '20c) - xp sin Sp cos(0i; - 26'^)] 

^ G"+ + R^bG'- — Rb [vd cos 5^ cos(0b — 26'c) + x/j cos 6b sin(0B — 2^c)] ' 

while 

. Rb [yp sin 6b sin(0B - 29^) + xp sin 6d cos(0b - 26'e)] 

^ R\G'^ + G'_ — -Rb [i/d cos (5b cos(0b — 2^c) — a;/) cos 6b sin((/)B — 2^c)] 

Here G"|_ _ = rC^^- = (1 + ^° ^^ ^ Vd ^d ^ respectively. The parameters xp and yp are 
subjected to stringent experimental bounds in case of ^c = 0,[21]: x\, + y]^ < (6.7 x 10~^)^. 
For non-vanishing dc-, this bound is no longer valid. However it is believed that in general 
Xp ^ yp ^ 10^^. In this case, it is clear that G'_,_ ~ 1 and G'_ ~ 10~^ which imply that 

Acp, - 10-2 xRb- 0(10-=^), (91) 

and 

AcP.^^^^OiQ.l). (92) 

Kb 

From these results, it is remarkable that the effect of D^ — Z)° mixing breaks the usual 

relation between the CP asymmetries Acp^ = A^ and AcP2 = Aj^: A+ ^ — A_ which 

is satisfied in the SM and SUSY models with negligible D^ — D^ mixing, as we have 

emphasized in the previous sections. As an example to show how natural to obtain in this 

case CP asymmetries of order their central vales of the experimental results in Eq.(l), let 
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us consider Rb ~ 0.15, x/) ~ 3 x 10 ^, y ~ 5 x 10 ^, 5^ ~ vr and 0^ ~ 6'c — 7r/4. One 
case easily find that 



i-CPi 



0.002, Acp^ ~ 0.3, (93) 



It is interesting to note that these values of the CP asymmetries depend on the CP 
violating SM phase 7 and the SUSY phase in the b ^ u transition 0i, which contribute 
together to (pB as in Eq.(69), in addition to the V^ — Z)° mixing phase 6c- Therefore, the 
determination of the angle 7 relies on the new SUSY phases 0i and Oc- This confirms the 
fact the our determination of the SM angel might be infiuenced by a new physics effect. 

6 Conclusions 

In this paper we have studied supersymmetric contributions to B^ -^ D^K^ and B^ — > 
D^K~ processes. We have shown that in the SM, the branching ratios Rcp± of these 
processes are within the experimental range. However the CP asymmetry Acpj^ is below 
its Icr experimental lower bound and the value of Acp_ is always negative. We have 
performed a model independent analysis of the gluino and chargino contributions tob ^ u 
and b ^ c transitions. We have used the mass insertion approximation method to provide 
analytical expressions for all the relevant Wilson coefficients. 

The D^ — D^ mixing experimental limits imply strong constraints on the mass inser- 
tions {S\b)i2 which affect the dominant gluino contribution to B^ -^ DK^ . We have 
revised these constraints and took them into account. We showed that in case of negligible 
D^ — Z)° mixing, it is possible to overcome these constraint and enhance the SUSY results 
for the CP asymmetries in B^ -^ DK~ if one consider simultaneous contributions from 
more than one mass insertion. In this case, the Acp^ becomes within la experimental 
range. However, with a large D^ — D^ mixing, one finds a significant deviation between 
the two asymmetries and it becomes natural to have them of order the central values of 
their experimental measurements. 

In general. We have emphasized that SUSY CP violating phases may contribute sig- 
nificantly to the CP asymmetries in B~ -^ DK~ and therefore, they may affect our 
determination for the angle 7 in the unitary triangle of the CKM mixing matrix. 
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